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ABSTRACT 


Three different multiscale models, based on the method of cells (generalized and 
high fidelity) micromechanics models were developed and used to predict the elastic 
properties of C/C-SiC composites. In particular, the following multiscale model- 
ing strategies were employed: Concurrent multiscale modeling of all phases using 
the generalized method of cells, synergistic (two-way coupling in space) multiscale 
modeling with the generalized method of cells, and hierarchical (one-way coupling 
in space) multiscale modeling with the high fidelity generalized method of cells. The 
three models are validated against data from a hierarchical multiscale finite element 
model in the literature for a repeating unit cell of C/C-SiC. Furthermore, the multi- 
scale models are used in conjunction with classical lamination theory to predict the 
stiffness of C/C-SiC plates manufactured via a wet filament winding and liquid sili- 
con infiltration process recently developed by the German Aerospace Institute. 


INTRODUCTION 


Advanced ceramic matrix composites (CMCs) have a wide range of applicabil- 
ity due to the extremely high operating temperatures (> 1400° C) the materials can 
sustain, as compared to polymer matrix composite (PMC) counterparts, and the im- 
proved toughness of CMCs over monolithic ceramics. These attractive features, along 
with the high specific stiffness and strength, make CMCs ideal for various applica- 
tions. Moreover, the material behaves more like a structure than a typical material 
because of the numerous phases, and as such, there is flexibility in designing the ma- 
terial itself. In addition, residual stresses and/or processing induced cracks and voids 
usually develop during the manufacturing which affect the performance of the CMC. 
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These manufacturing induced phenomena can be controlled through processing, fur- 
ther expanding the design space for the CMC, and thus, enabling the ability to tailor 
the material for specific performance criteria, as warranted by a structural design. 

A relatively new CMC system of particular interest consists of tows containing 
carbon (C) fibers embedded in amorphous C and surrounded by a silicon carbide 
(SiC) matrix (C/C-SiC) [1, 2]. C/C-SiC has found utility as an emerging technology 
for a wide range of applications for various reasons. The ability to tailor and function- 
ally grade the microstructure through processing and manufacturing enables a wider 
design space of C/C-SiC structures, and improves the ability to meet structural per- 
formance specification [3, 4]. Moreover, C/C-SiC systems exhibit desirable elastic, 
thermal, fracture, and damage tolerance properties [5, 6, 7]. Finally, the production 
of C/C-SiC is relatively inexpensive, particularly compared to other CMC systems 
[2, 8, 9]. 

Due to the hierarchical geometric structure and the direct link of the composite 
performance to the microstructure, micromechanics or multiscale modeling must be 
employed to accurately capture the response of C/C-SiC structures. Predictive mi- 
cromechanics and multiscale models of the thermomechanical behavior of C/C-SiC 
can be effective tools for taking full advantage of the customizable nature of the C/C- 
SiC material system. Blacklock and Hayhurst [10] developed a hierarchical multi- 
scale framework (microscale model is homogenized a priori and output is utilized as 
input at the meso- or macroscale) for predicting the elastic properties of C/C-SiC uti- 
lizing the finite element method (FEM). Furthermore, several authors have developed 
analytical micromechanical models to predict the non-linear mechanical response due 
to damage evolution of C/C-SiC [11, 12, 13, 14, 15, 16]. 

It has been demonstrated that filament winding is a promising manufacturing 
technique, as compared to the use of textile preforms, to attain reinforcement in ar- 
bitrary directions within CMC structures [9, 17]. A novel manufacturing method 
utilizing wet filament wind unidirectional carbon fiber perform and liquid silicon in- 
filtration (LSI) has been develop by the German Aerospace Institute (DLR) to manu- 
facture state-of-the-art C/C-SiC structures [18]. The filament winding manufacturing 
process offers increased flexibility in two ways. First, the lay-up of the fiber pre- 
form can be tailored to the meet local thermomechanical requirements of the struc- 
ture. Second, the C/C-SiC microstructure is a function of the filament winding angle 
[18]. Recently, filament wound C/C-SiC has been used to successfully manufacture 
a small-scale rocket engine nozzle extension [19]. 

The analytical micromechanics models lack the ability to directly capture dis- 
order, or non-uniformity, in the arrangement of the constituents of the composites, 
and often cannot resolve gradients in the local fields (stress and strain), properties, 
and/or material state. These deficiencies can be overcome with fully-numerical FEM 
or DEM. Yet, the computational intensity of these models impede their applicability 
within a synergistic (hierarchical in space, but concurrent in time) multiscale frame- 
work. In addition, little (if any) multiscale modeling of filament wound C/C-SiC 
plates or structures is present in the current literature. Thus, the goal of this cur- 
rent research effort was to develop and evaluate several multiscale models of fila- 
ment wound C/C-SiC which utilize semi-analytical micromechanics theories that are 
amenable to a synergistic implementation. 


The method of cells (MoC) family of semi-analytical micromechanics theories, 
including the generalized method of cells (GMC) and high fidelity generalized method 
of cells (HFGMC), offer greater fidelity (including local gradients in the fields, prop- 
erties and material state) than analytical micromechanics models and are more effi- 
cient than numerical models [20]. Thus, the MoC theories are suitable for imple- 
mentation within a synergistic multiscale modeling scheme. Three different multi- 
scale modeling strategies were employed: 1. Concurrent multiscale modeling of all 
phases using GMC, herein referred to as “GMC,” 2. Synergistic (two-way coupling 
in space) multiscale modeling with GMC, referred to as “MsGMC” [21], 3. Hier- 
archical (one-way coupling in space) multiscale modeling with HFGMC, denoted 
as “MsHFGMC”. The three models are validated against results from a hierarchical 
multiscale FEM model for a repeating unit cell of C/C-SiC in the available literature 
[10]. Furthermore, the multiscale models are used in conjunction with classical lam- 
ination theory to predict the stiffness of C/C-SiC plates manufactured via a filament 
winding process recently developed by DLR [22]. 

The organization of this paper is as follows. The various novel multiscale frame- 
works developed for C/C-SiC are presented along with the underlying micromechan- 
ics theories. Then, validation results of the multiscale models as compared to data 
from existing computational results [10] are given. Finally, the multiscale models 
are used to predict the elastic properties of filament wound plates and compared to 
experiment. 
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Figure 1. Micrograph of C/C-SiC displaying RUC with distinct phases: C fibers within 
amorphous C matrix (C/C bundle), SiC matrix, and free Si[22]. 
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Figure 2. Parameterization of C/C-SiC RUC. 


MULTISCALE MODELING OF C/C-SIC RUC 
Parameterization of C/C-SiC RUC 


C/C-SiC composites contain three dominant phases (C fiber, amorphous C ma- 
trix, and SiC matrix), excluding free silicon (Si), existing at various length scales, see 
Figure | [22]. In addition, a repeating unit cell (RUC) within the C/C-SiC microstruc- 
ture can be readily identified, marked with a green box in Figure 1. These features 
facilitate the modeling of C/C-SiC via a multiscale strategy. In order to represent 
the microstructure of the filament wound C/C-SiC composite with an RUC, key fea- 
tures of the microstructure must be identified, idealized, and parameterized. Figure 2 
displays the RUC obtained from Figure 1. The C fibers are shown in dark blue, the 
amorphous C matrix is shown in green and the SiC matrix is shown in yellow. The 
combination of the C fibers and amorphous C matrix is referred to as a “C/C bundle.” 
The key parameters within the C/C-SiC RUC are described in Table I. Once these 
parameters are quantified, along with the volume fraction of each phase, the geom- 
etry of the C/C-SiC RUC is completely defined. Following Blacklock and Hayhurst 
[10], the geometric parameters for the C/C-SiC RUC are given in Table I, and the 
volume fractions of the constituents are given in Table II. Through homogenization 
and application of periodic boundary conditions, the response of the RUC represents 
that of an infinitesimal point in a continuum at a higher scale. Therefore, there is no 
dependence on length scale, only relative volume fractions of the constituents within 
the RUC. Hence, the geometric parameters given in Table I are dimensionless. 


Concurrent Multiscale Model of C/C-SiC RUC with GMC 


The doubly-periodic formulation of GMC was used to calculate the effective elas- 
tic properties of a C/C-SiC RUC, wherein all present phases are modeled concurrently 


TABLE I. LIST AND DESCRIPTION OF KEY MICROSTRUCTURAL PARAMETERS 

USED TO DEFINE C/C-SIC RUC. 
Parameter Description Value 
Abundle X_ dimension of C/C bundle 1.0 
oinais x3 dimension of C/C bundle 0.15 
dy Diameter of C fiber 0.029 
Of Fiber spacing 0.0012 
tsic Thickness of SiC matrix 0.033 
€2 Distance from x2 boundary to nearest fiber 0.0011 
€3 Distance from x3 boundary to nearest fiber 0.0007 
Ny Number of fibers in x» direction 27 
n3 Number of fibers in x3 direction 4 


TABLE II. LIST AND DESCRIPTION OF CONSTITUENT VOLUME FRACTIONS 
USED TO DEFINE C/C-SIC RUC. 


Parameter | Description Value 
vs, C fiber volume fraction 0.49 
ug Amorphous C matrix volume fraction 0.16 
Ugic SiC matrix volume fraction O35 
vc/C Volume fraction of C/C bundle 0.65 
ue, CO Volume fraction of C fiber in C/C bundle 0.76 
VEC Volume fraction of amorphous C matrix in C/C bundle | 0.24 


at a single scale. GMC, first developed by Paley and Aboudi [23], offers a very com- 
putational efficient means of calculating the local fields within and effective proper- 
ties of a composite RUC. A brief summary of the GMC fundamentals is presented in 
this section; the reader is referred to Aboudi et al. [20] for complete details on GMC. 
With GMC, an RUC must be identified within the composite microstructure, as 
described in the previous section for the C/C-SiC composites. The RUC is then dis- 
cretized into a number of rectangular subcells. The C fiber reinforcement in the C/C 
bundles are continuous in the axial, x,-direction. Thus, it is assumed that the RUC 
extends infinitely in that direction, and the RUC is subjected to a state of general- 
ized plane strain. Therefore, the doubly-periodic formulation of GMC is presented 
herein. Subcell indices 6 and ¥ are used the location of the subcell within the RUC 
in the local x2-, and x3-directions, respectively. Each subcell can be occupied by a 
single constituent (obeying the constitutive law of that constituent), but there is no 
limit on the number of different constituents that may be present in the RUC. The 
displacement fields u in each subcell {y are approximated to be linear. 


uw27) = wi) 4 gL) pv) zfs EM pO (1) 


where w'*”) are the displacement components of the center of the subcells, and pe” 
and ap? ” characterize the linear dependence of the displacements on the local subcell 
coordinates %2 and 73. Displacement and traction continuity is enforced in an aver- 


age, or integral sense at each of the subcell interfaces and the periodic boundaries 
of the RUC. These continuity conditions are used to formulate a strain concentra- 
tion matrix A‘), which gives all the local subcell strains (e(87)) in terms of the six, 
global, average, applied strains €. 


6h) = ALINE (2) 


Noteworthy, even though the formulation is doubly-periodic, the global RUC, and 
local subcell, stress and strain fields are fully-3D. 

Once the strain concentration matrix is formulated and the local subcell strains 
have been calculated, the local subcell stresses (0 “°”) can be calculated using the 
local constitutive law and the local subcell strains. 


o FY) = Clb) (1) (3) 
where C7) is the local subcell stiffness, €{°”) contains the six local subcell strains. 
Eq. (3) can be adjusted to account for both thermal and inelastic strains. However, 
the analyses considered here are limited to linear elastic materials. Volume weighted 
averaging can be employed to calculate the average global stresses o@ in terms of the 


local subcell stresses. 
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Finally, the global, homogenized RUC stiffness C* can be calculated. 
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where Ng, and NV, are the total number of subcells in the x2-, and x3-directions, 
respectively, H, and L are the RUC dimensions in the x2-, and x3-directions, hg, and 
L, are the local subcell dimensions in the x2-, and x3-directions for subcell 3. 

The RUC used to represent the C/C-SiC microstructure, in a concurrent manner, 
is shown in Figure 3a, and the associated subcell grid is given in Figure 3b. C fiber 
subcells are represented in dark blue, amorphous C matrix in green, and SiC matrix 
in yellow. This color scheme is used throughout this paper. While the arrangement 
of fibers in the physical composite is random, a square-packed architecture was cho- 
sen as an idealization. Modeling the C/C-SiC with a random microstructure would 
require statistical considerations when calculating the effective properties. There re- 
mained some flexibility in other RUC parameters including number of fibers, fiber 
diameter, and fiber spacing. A thorough parametric study was conducted to arrive at 
the RUC displayed in Figure 3. A uniform distribution of fibers is desired because 
the lack of coupling between normal and shear stresses, inherent in the formulation of 
GMC, leads to error in the calculation of the local stresses and homogenized proper- 
ties [24]. Minimization of the difference between e2 and e3, see Figure 2, represents 
the most uniform distribution of fibers. It was deemed that the most uniform config- 
uration (smallest difference between ez and e3) could be achieved with four fibers in 
the x3-direction (n3) and 27 fibers in the x-direction (n2), as shown in Figure 3. 


(a) No subcell grid. 


(b) Subcell grid displayed. 


Figure 3. Concurrent GMC model of C/C-SiC RUC with C fiber (dark blue), amorphous C 
matrix (green), and SiC matrix (yellow). 


Synergistic Multiscale Model of C/C-SiC RUC with MsGMC 


A recursive formulation of GMC, MsGMC, is utilized to model the C/C-SiC RUC 
across multiple length scales in a synergistic manner (i.e., hierarchical in space but 
concurrent in time)[21]. MsGMC can admit any general number of scales, or levels, 
denoted by 7. Thus, the strain concentration, local stress, and homogenized stiffness 
relationships given in Eqs. (2)-(5) hold for each level 2 in the analysis. 
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The separate scales are linked by equating the homogenized strains, stresses, and 
stiffness at Level 7 (€;, &;, C*, respectively) to the local strains, stresses, and stiffness 
within a given subcell at the next highest scale, Level i-1 (677, 08), C(8, 
respectively). 


a= ef) (10) 


(a) Level 1. Mesoscale: SiC matrix (yellow) and homogenized C/C bundle (gold). 


(b) Level 2. Microscale: C fiber 
(dark blue) and amorphous C matrix 
(green). 


Figure 4. MsGMC RUCs and associated subcell grids. 
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Transformations can be performed such that the lower scale RUC is oriented with 
respect to a different coordinate system than that of the subcell it is linked to at the 
higher scale [21]. However, the transformations were omitted from the formulation 
presented here because all scales utilize the same coordinate-basis vectors. Eqs. (10) 
- (12) can be used in conjunction with Eqs. (6), (8), and (9) to obtain the homogenized 
strains, stresses, and stiffness at any scale. Localization and homogenization across 
any general number of scales is automated. 

The C/C-SiC composite is modeled using two levels within MsGMC and the 
hierarchy shown in Figure 4. At the mesoscale (Level 1) displayed in Figure 4a, the 
RUC for the C/C-SiC composite is composed of SiC matrix, shown in yellow, and the 
homogenized C/C bundle, shown in gold. Each subcell representing the C/C bundle 
is linked to a microscale (Level 2) RUC, presented in Figure 4b, containing a single 
C fiber within amorphous C matrix depicting a square packed architecture. 

The geometry used to represent the fiber, shown in Figure 4b, was chosen after 
a thorough convergence study. For a fixed geometry, the local and global results 
are completely insensitive to refinements in the subcell grid used to represent the 
geometry because of the lack of normal-shear coupling in GMC. Thus, the grids 


shown in Figure 4 represent the minimum level of discretization needed to represent 
the chosen geometry. 


Hierarchical Multiscale Model of C/C-SiC RUC with MSHFGMC 


GMC offers a computationally efficient platform for homogenization of compos- 
ite microstructures, and provides reasonable accurate results for effective properties 
and other averaged quantities. Conversely, HFGMC offers more accurate local fields 
at a computational expense, but attains normal-shear coupling. A third model was 
devised, MSHFGMC, which utilizes the HFGMC micromechanics theory across mul- 
tiple scales to model the C/C-SiC composite . However, unlike MsGMC, the scale 
bridging has not yet been automated (synergistic). Instead, homogenization is per- 
formed in a “bottom-up,” one-way coupling hierarchical scheme. It should be noted 
that, there exists no difference between the homogenized composite stiffness results 
obtained from hierarchical and synergistic multiscale models for a linear analysis, 
like that performed herein. In this section, a brief overview of the HFGMC microme- 
chanics theory is given, and the reader is referred to Aboudi, et al. for more details 
[20]. 

HFGMC utilizes quadratic approximations of the displacement fields in each sub- 
cells at Level 7. 
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averaged displacements, and wee a (when both m # 0 and n ¥ 0) characterize the 
quadratic dependence of the Sieslacements on the local subcell coordinates. The un- 
known terms we are determined through satisfaction of equilibrium, enforcement 
of periodic boundary conditions and continuity of displacements and tractions at the 
subcell interfaces. As with GMC, continuity conditions are satisfied in an average 
(integral) sense. 

Imposing these conditions at each level (scale) yields a linear set of algebraic 
equations. 


where €; are the externally applied average strain components, w! are the volume- 


K,U; =f, (14) 


The K; matrix contains geometric and stiffness terms for the individual subcells at 
Level 2, the displacement vector U; is composed of the unknown displacement coef- 
ficients (i.e., components of W, 6 os) at Level 7, and the force vector f; consists of the 
applied average strains €; ay eel i. A localization relationship, similar to Eq. (6) 
for GMC, between the applied average strains and the local subcell strains at Level 7 
can be obtained by solving Eq. (14). 

eR) = APE, (15) 
where Ae is the HFGMC strain concentration matrix at Level 7. Once the strain 
concentration matrix has been calculated, the homogenized stiffness at Level 7 can be 


readily determined. 
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Since the scale bridging is not automated, the local fields in MsHFGMC are not 
calculated, and handshaking across the scales is achieved in a bottom-up hierarchical 
manner with Eq. (12). 

The MsHFGMC model operates at the same two length scales, defined in Figure 
5, as the MsGMC model described previously. The same geometry at Level | was 
used to represent the SiC matrix and homogenized C/C bundle in both MsGMC and 
MsHFGMC models. Unlike GMC, HFGMC is sensitive to refinements in the subcell 
grid. A convergence study was performed to determine the appropriate subcell grid 
density for Level 1 shown in Figure 5a. Finally, a study on both refinement in the fiber 
geometry and subcell gird were performed to arrive at the RUC for the C/C bundle 
displayed in Figure 5b. 


NUMERICAL VALIDATION OF MULTISCALE MODELS OF C/C-SIC RUC 


In this section predictions for the elastic properties of C/C-SiC obtained using 
the three models presented in the preceding section are validated against results from 


(a) Level 1. Mesoscale: SiC matrix (yellow) and homogenized C/C bundle (gold). 


(b) Level 2. Microscale: C fiber 
(dark blue) and amorphous C matrix 
(green). 


Figure 5. MsHFGMC RUCs and associated subcell grids. 


TABLE III. ELASTIC PROPERTIES FOR C FIBER AND AMORPHOUS C MATRIX 
USED IN MULTISCALE MODELS. 


C Fiber Value | Amorphous C Matrix | Value 
Evi (GPa) | 294.0 | EG» (GPa) 62.6 
Ess (GPa) | 18.4 | vm 0.3 
vi 0.2 | Gon (GPa) 24.0 
Usd 0.25 
Gi (GPa) | 8.96 


TABLE IV. ELASTIC PROPERTIES FOR SIC MATRIX USED IN MULTISCALE MOD- 
ELS. 


Damaged SiC Matrix | Value | Undamaged SiC Matrix | Value 
Bei (MPa) 0.001 | EC (GPa) 15.0 
EXiC™ (GPa) 15.0 | pStCm 0.3 
ea 2.0E-8 | G°*Cm (GPa) 577 
Dada 0.3 
GtCm (GPa) 577 


a bottom-up, hierarchical multiscale FEM model of C/C-SiC tows within a woven 
composite [10]. The RUCs modeled in this work utilized the same geometrical and 
volume fraction data (given in Tables I and II), although the intent of the current work 
is to predict the elastic properties of C/C-SiC plates manufactured using a wet fila- 
ment winding technique. In addition, the elastic properties for the C fiber, amorphous 
C matrix, and SiC matrix constituents were also taken from Blacklock and Hayhurst 
[10] and given in Tables III and IV. 

In Blacklock and Hayhurst [10] it was stated that cracks developed in the SiC 
matrix during processing. As such, the value for the SiC stiffness in the direction 
aligned with the fiber axis £7/©™, used in the FEM model was chosen to be zero . 
There is no mention of how the other properties were treated; so, it is assumed that 
these other properties did not change from the undamaged isotropic values, yielding 
an anisotropic stiffness matrix for the cracked SiC. However, in the doubly-periodic 
formulation of GMC and HFGMC, used for the analyses presented here, it is assumed 
that the components of the stiffness matrix are transversely isotropic. Thus, if the 
value for the axial Poisson’s ratio v7}°" was not set to a very low value, as listed in 
Table IV, some of the resulting Young’s moduli for the homogenized C/C-SiC RUC 
would be non-physical (i.e., negative). 

The results for all six engineering constants for the C/C-SiC RUC predicted using 
GMC, MsGMC, and MsHFGMC are presented in Table V along with the percentage 
difference from the FEM predictions of Blacklock and Hayhurst [10]. MsHFGMC 
predicted elastic properties closest to the FEM analysis for nearly all properties, fol- 
lowed by MsGMC and the GMC. The predictions from the models matched most 
closely to the FEM for the Young’s moduli, exhibiting a ~0-10% difference, and the 
predictions for the shear moduli were within ~3-18% of the FEM. A significant dis- 


TABLE V. ELASTIC PROPERTIES OF C/C-SIC RUC PREDICTED WITH MULTISCALE 
MOC MODELS, INCORPORATING DAMAGED SIC, COMPARED TO MULTISCALE 
FEM CALCULATIONS FROM BLACKLOCK AND HAYURST, 2012. 
Property | GMC | % Difference | MsGMC | % Difference | MsHFGMC | % Difference 
from FEM from FEM from FEM 
Ey, (GPa) | 155.2 0.6 155.2 0.6 155.2 0.6 
E92 (GPa) | 20.5 6.4 20.7 as) 213 Zeh 
E33 (GPa) | 19.6 9.3 20.1 7.0 21.6 4.9 
Vi2 0.167 34.4 0.168 34.0 0.173 32.0 
V3 0.143 45.2 0.146 44.2 0.149 43.1 
V23 0.272 12.4 0.268 13.8 0.272 25 
G23 (GPa) | 7.45 8.9 7.45 8.9 7.94 2.8 
G13 (GPa) | 8.34 18.3 8.44 17.2 8.60 15.7 
Gig (GPa) 8.96 1221 9.01 11.6 275 9.3 
TABLE VI. ELASTIC PROPERTIES OF C/C-SIC RUC PREDICTED WITH MULITI- 
SCALE MOC MODELS, INCORPORATING UNDAMAGED SIC, COMPARED TO MUL- 


TISCALE FEM CALCULATIONS FROM BLACKLOCK AND HAYURST, 2012. 

Property | GMC | % Difference | MsGMC | % Difference | MsHFGMC | % Difference 

from FEM from FEM from FEM 
Ey, (GPa) | 160.0 3.9 160.0 3.9 160.0 3.8 
E99 (GPa) | 21.0 4.0 74 3.0 21:9 0.3 
E33 (GPa) | 20.5 5.0 21.1 2.6 21.6 0.2 
V2 0.251 1.3 0.251 1.2 0.254 0.0 
V3 0.255 22 0.258 1.2 0.261 0.0 
V3 0.312 0.3 0.308 1.1 0.312 0.2 
G3 (GPa) | 7.45 8.87 7.45 8.9 7.94 2.8 
G13 (GPa) | 8.34 18.3 8.44 172 8.60 15.7 
G'j2 (GPa) | 8.96 1221 9.01 11.6 9.25 93 


crepancy (~12-45%) between the Poisson’s ratios predicted here and that from the 
FEM model was observed. This can be attributed to the value for axial Poisson’s 
ratio for the damaged SiC matrix v>}“™ that was needed, see Table IV, because of the 
assumption of transverse isotropy in the double-periodic MoC-based formulations. If 
it is assumed the SiC matrix in the MoC-based models is undamaged and isotropic 
properties are used, a drastic reduction to 0-2% in difference between the Poisson’s 
ratios of the C/C-SiC RUC predicted using the current models and the FEM model in 
the literature is observed, as shown in Table VI. 

All multiscale simulations were executed on a single 2.40 GHz Intel Xeon E5- 
2630 v3 processor (cpu). Figure 6 shows the analysis time versus total number of 
subcells for each model. The run time increases with the total number of subcells in 
the model with MsGMC being the fastest. However, all analyses completed in less 
than 10 seconds. 


Run Time (Seconds) 


0 2000 4000 6000 8000 
Total Number of Subcells 


Figure 6. Analysis time versus total number of subcells for three multiscale models of C/C- 
SIC. 


The results given in Tables V and VI show that GMC, MsGMC, and MsHFGMC, 
can predict elastic properties for a C/C-SiC that are comparable to previously pub- 
lished FEM results and serve as a numerical validation of these methods. Herein, 
subsequent results will be normalized to experimental data or baseline results from 
the MoC-based models. 


RESULTS AND DISCUSSION 


A focus of this paper is to utilize the MoC-based models (GMC, MsGMC, and 
MsHFGMC) to predict the elastic properties of C/C-SiC plates, manufactured by 
Breede and colleagues using a wet filament winding and LSI technique, where it has 
been shown that the C/C-SiC microstructure is dependent on the filament winding 
angle [18]. Unfortunately, significant data on the geometry and volume fraction of 
the constituents within the filament wound C/C-SiC, needed to create a representa- 
tive RUC, was not available. Therefore, the same RUC used for the predictions in the 
previous section with the geometry and constituent volume fractions given in Tables I 
and II for a tow within a woven C/C-SiC material are used to predict the elastic prop- 
erties of filament wound C/C-SiC. The elastic properties for the constituent materials 
presented in Tables II and IV are also used. To model the filament wound C/C-SiC 
plates, the homogenized response the C/C-SiC RUC is used to represent the behavior 
of the layers within a [6/-0] laminate using classical lamination theory (CLT), where 
6 is the filament winding angle [25]. The available reports did not comment if the 
filament wound plates were manufactured to have a symmetric configuration. Yet, it 
is assumed here that the plates are symmetric to avoid coupling between axial normal 
loads and shear loads or moments during coupon testing. 

Three different filament winding configurations ([0°/90°], [15°/-15°], and [45°/- 
45°]) were manufactured and tested by Breede et al. [22]. Figure 7 displays the 
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Figure 7. Prediction of Young’s modulus F,,.,, of filament wound C/C-SiC at various winding 
angles compared to experiment from Breede et al. [22]. 


predicted values for the Young’s moduli in the loading direction F’,,, of the coupons, 
normalized to the experimental mean, for three different filament winding angles 
(#). The stiffness predictions for the [0°/90°] and [45°/-45°] configurations showed 
excellent agreement (< 3% error) with the experimental data, and are fairly consistent 
across all three multiscale models. Conversely, all the models exhibit ~20% error 
when predicting the stiffness of the [15°/-15°] plate. Substantial differences between 
the microstructures of the [45°/-45°] and [15°/-15°] plates were reported by Breede 
et al. [22]. Thus, it is inferred that the idealization of the microstructure of the tows 
within the woven C/C-SiC more closely represents that of the [0°/90°] and [45°/-45°] 
filament wound C/C-SiC, whereas new geometric information would be required to 
accurately predict the properties of the [15°/-15°] filament wound C/C-SiC using the 
MoC-based multiscale models. 

It is expected that multi-axial stress states will be present throughout the filament 
wound structure during service. Thus, the transverse stiffness E,,, and shear stiff- 
ness G,,,, of the filament wound plates are also important. Figure 8 shows predictions 
of E,,, Ey,, and G,, for all three filament winding angles obtained with all three 
multiscale models. From these predictions, it can be deduced that the [0°/90°] has 
an intermediate axial and transverse stiffness but low shear stiffness, as compared to 
the other lay-ups. The [15°/-15°] exhibits the highest axial stiffness, but low trans- 
verse and shear stiffness. Finally, the [45°/-45°] offers the highest shear stiffness. 
These stiffness predictions can be utilized in the preliminary design phase of a C/C- 
SiC structure to meet the stiffness requirements throughout the structure under local, 
multi-axial stress states. In addition, with the appropriate data, the microstructure in 
the different lay-ups can also be modeled as a function of the filament winding angle. 
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Figure 8. Prediction of elastic properties of filament wound C/C-SiC at various winding 
angles 


CONCLUSIONS 


Three models were developed for predicting the elastic properties of C/C-SiC 
manufactured using wet filament winding and LSI. These models utilized the GMC 
and HFGMC micromechanics theories across multiple length scales. The first model 
consisted of a concurrent length scale representation of all three phases within the 
C/C-SiC composite RUC an utilized GMC for localization and homogenization. The 
second scale incorporated two separate length scales: a mesoscale containing the 
SiC matrix and homogenized C/C bundle, and a microscale containing the C fiber 
and amorphous C matrix. Communication between the scales was accomplished in 
a two-way coupling, synergistic manner (i.e., concurrent in time, but hierarchical in 
space) through localization and homogenization with GMC. The final model utilized 
HFGMC to calculated the homogenized properties of the microscale RUC for use in 
the constitutive relationship of the C/C bundle at the mesoscale via a bottom-up, one 
way coupling, hierarchical multiscale framework. 

All three multiscale models were used to predict the effective properties of an 
RUC within a woven C/C-SiC composite, and showed excellent agreement with anal- 
ogous FEM predictions available in the literature. Furthermore, the run time for all 
simulations was less than 10 seconds using a single CPU. The same material proper- 
ties and geometric parameters for the RUCs were used to predict the elastic properties 
of filament wound C/C-SiC with three different winding angles. Predictions, by all 
three models, for two out of the three configurations exhibited excellent agreement 
with experimental data. The error exhibited in the predictions for the third configura- 
tion [15°/-15°] support the micro-graphical evidence that the C/C-SiC microstructure 
changes depending on the filament winding angle. To capture this the RUCs used for 
this particular angle would need to be altered based upon measurements taken at the 
microscale. 

The effect of both free Si and C matrix porosity on the elastic properties of a 
C/C-SiC composite can be resolved by adjusting the effective properties of the ho- 
mogenized representation of SiC and amorphous C. However, this will not produce 
the same effect on the local fields within the composite. The local fields are the main 
driving force for any damage initiation or inelasticity that may arise and must be 
considered if accurate non-linear predictions are required. The multiscale analysis 
presented in this work shows that the addition of free Si or C matrix porosity tends 
to reduce the average effective stress within the C matrix phase of the composite. 
While the free Si also reduces the maximum effective stress, locally within the C ma- 
trix, the presence of pores increases this value and serves as a stress riser. Full-field 
micromechanics theories are necessary to capture such effects. 

It has been demonstrated that, multiscale modeling frameworks utilizing GMC 
and HFGMC are in place for predicting the elastic properties of filament wound C/C- 
SiC at every length scale. Additional phases can be included in the models, and it is 
possible to represent the microstructure as a function of the filament winding angle 
(although this was not demonstrated here). With this framework in place, full-scale, 
filament wound C/C-SiC structures can be tailored and optimized to meet stiffness 
and strength requirements in a computationally efficient manner. In addition, other 
aspects of the filament wound C/C-SiC composite can be modeled including thermal 


properties and effects, processing, and damage evolution. 
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